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Abstract 

The random greedy algorithm for constructing a large partial Steiner-Triple-System is defined 
as follows. Begin with a complete graph on n vertices and proceed to remove the edges of 
triangles one at a time, where each triangle removed is chosen uniformly at random out of all 
remaining triangles. This stochastic process terminates once it arrives at a triangle-free graph, 
and a longstanding open problem is to estimate the final number of edges, or equivalently the 
time it takes the process to conclude. The intuition that the edge distribution is roughly uniform 
at all times led to a folklore conjecture that the final number of edges is v?'/'^+°^^'> with high 
probability, whereas the best known upper bound is -nJ / It is no coincidence that various 
methods break precisely at the exponent 7/4 as it corresponds to the inherent barrier where 
co-degrees become comparable to the variations in their values that arose earlier in the process. 

In this work we significantly improve upon the previous bounds by establishing that w.h.p. 
the number of edges in the final graph is at most n^~^^^°'^^\ Our approach relies on a system 
of martingales used to control key graph parameters, where the crucial new idea is to harness 
the self-correcting nature of the process in order to control these parameters well beyond the 
point where their early variation matches the order of their expectation. 



1 Introduction 

We consider the random greedy algorithm for triangle-packing. This stochastic graph process begins 
with the graph G(0), set to be the complete graph on vertex set [n], then proceeds to repeatedly 
remove the edges of randomly chosen triangles (i.e. copies of K^) from the graph. Namely, letting 
G{i) denote the graph that remains after i triangles have been removed, the {i + l)-th triangle 
removed is chosen uniformly at random from the set of all triangles in G{i). The process terminates 
at a triangle-free graph G{M). In this work we study the random variable M, the number of 
triangles removed before obtaining a triangle-free graph (or equivalently, the number of edges in 
the final triangle- free graph, which is (2) — 3M). 
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This process and its variations play an important role in the history of combinatorics. Rodl |13j 
proved the Erdos-Hanani conjecture — which posits the existence of large partial Steiner systems, 
collections of t-seis with the property that no k-set is a subset of more than one set in the collection 
— in the early 1980's by way of a randomized construction that is now known as the Rodl nibble. 
This construction is a semi-random variation on the random greedy packing process defined above. 
It is semi-random in the sense that the desired object is constructed in a sequence of substantial 
pieces, where the proof of the existence of each piece is an application of the probabilistic method. 
Such semi-random constructions have been successfully applied to establish various key results in 
combinatorics (see pQ for an early application of this approach and [3] and [T^] for further details). 
We note in passing that semi-random techniques have been refined to show the existence of partial 
Steiner systems that are nearly as large as allowed by the simple volume bound, see [11] and |17j . 
In particular, Alon, Kim and Spencer [2] used such techniques to prove the existence of a collection 
of edge-disjoint triangles on n vertices that covers all but 0(n^/^ log"^^^ n) of the (2) edges of the 
complete graph. 

Despite the success of the Rodl nibble, the limiting behavior of the random greedy packing 
process itself remains unknown, even in the special case of triangle packing considered here. Recall 
that G{i) is the graph remaining after i triangles have been removed. Let E(i) be the edge set 
of G{i). Note that |£'(i)| = (2) — 3i and that E{M) is the number of edges in the triangle-free 
graph produced by the process. It is widely believed that the graph produced by the random 
greedy triangle-packing process behaves like the Erdos-Renyi random graph with the same edge 
density, hence the process should end once its number of remaining edges becomes comparable to 
the number of triangles in the corresponding Erdos-Renyi random graph (i.e., once |£^(M)| matches 
the order of (|£;(M)|/(2))^(3)). Here and throughout the paper we say that an event holds with 
high probability (w.h.p.) to denote that its probability tends to 1 as n — >• 00. 

Conjecture (Folklore). With high probability \E{M)\ = n^/'^+°('^\ 

Joel Spencer has offered $200 for a resolution of this question. It was shown by Spencer [15] in 
1995, and independently by Rodl and Thoma [l3] in 1996, that \E{M)\ = o{n^) w.h.p. Grable [TO] 
improved this bound to |£'(M)| < 

^11/6+0(1) ^-j^ 

adaptation of the Rodl nibble method and 
further sketched how similar arguments using more delicate calculations should extend this to a 
bound of rJ/^~^°W w.h.p. Wormald [18] later demonstrated how the differential equation method 
can also give nontrivial bounds for this problem (as well as generalizations of it), and namely that 
|£^(M)| < n^~57+°(^). Finally, in a companion paper [5] that introduced a differential equation 
approach to this process exploiting its self-correction nature, the foundations of the present work, 
the authors gave a short proof that |£'(M)| = 0{n'^^'^'^°^^^) w.h.p. 

It is important to note that the point at which there are roughly n^/^ remaining edges is a 
natural barrier in the analysis of this process. To illustrate this, notice that if the {i + l)-st 
triangle taken is abc then the change in the number of triangles in the graph once abc is removed 
is simply —\Na^b{i)\ — \Na,cii)\ — \^b,c{i)\ + 2, where Nu^v{i) denotes the common neighborhood of 
the vertices u,v £ [n] in the graph G{i). Hence, a natural prerequisite to analyzing this process 
is the understanding of the co-degrees \Nu^v\ for all u,v. Suppose for the sake of this discussion 
that early in the evolution of the process G{i) closely resembles the random graph with the same 
number of edges; that is, suppose G{i) is roughly the same as Gn,p where p = p{i) = 1 — 8^/(2) • If 
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this is the case when p is close to 1/2 (i.e. i is nearly (2) /6) then we expect the lA^^i^til's to be close 
to n/4 with variations as large as ^/n. If these variations in co-degrees persist to the point where 
p = n~^/^ (that is, i roughly (2) /3 — n''/^), where we expect the lA^^^^^I's themselves to be roughly 
n^/^, then these variations would be as large as their average value. Once this happens all control 
over co-degrees is lost, e.g. one could have all co-degrees with non-vanishing probability, or half 
of the co-degrees and the other half around n^/^, etc. In any case, if the variations in |A^M,t,| that 
develop early in the process are not somehow dealt with, one would expect the analysis to break 
down once n^/^ edges remain. Perhaps this is the reason that Wormald [18], who also treated this 
process with the differential equation method, stated that "some non-trivial modification would be 
required to equal or better Grable's result." 

In this work we exploit the self-correcting nature of the process in a system of carefully con- 
structed martingales which allows us to tighten the control over key graph properties over time and 
overcome the variations in their values that arise early in the process. Our main result is an upper 
bound on |i?(M)| that is significantly better than rJ^'^. 

Theorem 1. Consider the random greedy algorithm for triangle-packing on n vertices. Let M be 

the number of steps it takes the algorithm to terminate and let E[M) be the edges of the resulting 
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triangle-free graph. Then with high probability \E{M)\ = 0{n 2V2 log n). 

A key feature of our proof of Theorem [1] is an estimate for | A^^^^yl in which the variation decreases 
as the process evolves. We stress that estimates for random graph processes with this property are 
not commonly obtained by martingale arguments or the differential equation method. 

The remainder of the paper is organized as follows. In the next section we discuss our analysis 
of this process in more detail, listing the random variables that we track and the estimates on 
them that we are able to prove. The proof of our main result. Theorem 12.1^ follows in Section 3. 
Theorem [1] follows directly from Theorem 12.11 

2 Self-correcting Estimates 

Let {J-i) be the filtration given by the underlying process. We note in passing that our probability 
space is the set of all maximal sequences of edge-disjoint triangles on vertex set [n] with probability 
measure given by the uniform random choice at each step. 

For u,v,w S [n] define Nu = Nu{i) = {x £ [re] : xu G E{i)}, let A^^^^ = Ci and let 
Nu,v,w = n n Nu]. Our main interest is in tracking the number of triangles in G{i) and the 
variables 

Yu,v{i) = \NuA^)\ = \{x e W : e -^(^)> e • 
In the course of our argument we will also need to consider the variables 

Yu{i) = \Nu{i)\ and y«,„,^(i) = \Nu,v,w{i) \ ■ 

We begin by writing the one-step expected changes in our main variables of interest. For any 
random variable W let AVF be the one-step change AVF = W{i + 1) — W{i). Let Q{i) be the 
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number of triangles in G{i). We have 

E[Ay„,„ I Ti] 

IE[AQ I Ti] 



E 



Q 



xyzi^Q 



Q 



(2.1) 
(2.2) 



We use these one-step expected changes to relate the random variables to functions of a continuous 
'time' (following the approach to the differential equation method introduced in [1]). We choose 
the time-scaling t = t{i) = i/v? . Following the convention established in the Introduction we set 



p = p{i,n) = 1 



6i 



l-6t. 



n 



(2.3) 



Note that p can now be viewed as either a function of i or the continuous time t; we pass between 
these interpretations without further comment throughout the paper. Now, these choices yield the 
trajectories Yu^v{i) ~ 2/(0'^ ^'^d Q{i) ^ q(t)n where we set 



y{t)=p\t) 



and 



/(t) =/(t)/6. 



We can arrive at these equations by either deriving them from (j2.ip and (|2.2p and the assumption 
that the one-step changes in the trajectory are equal to the expected one-step change in the cor- 
responding random variable or by appealing to our Gn,p intuition. The companion paper [5] uses 
these variables alone to establish a bound of 0(n'^/^ log^''^ n) on the number of edges that survive 
to the conclusion of the algorithm. 

In order to achieve better precision, we introduce additional variables with the central goal of 
establishing an estimate for Y^^^ with variation that decreases as the process evolves. (For appli- 
cations of the differential equation method that exploit this kind of 'self-correcting' phenomenon, 
see [7] and |16j.) We would like to add random variables to our collection that will give us better 
control on the expression in the numerator of (|2.ip . the one-step expected change in Y^^y. To this 
end we take a closer look at this expression. We have 



xeNu,v 



(2.4) 



where 



p 



{xy G E{i) : x € Ny^y , y £ Nu\ Ny^y and y ^ v} 
{xy £ E{i) : x G Nu,y , y £ Ny\ Nu,y and y ^ u} 

E(i) n 



Ny 



We will use appropriate martingales to track each of these variables for every pair of vertices uv. 
Now, we expect these variables to be correlated with Yy,Yy and Yy^y. Indeed, we should have 



Ru,V ~ P {Yu Yy yyYy y , S y ^y ~ P {Yy Yy^y^Yy^.^ 



P 



Y 

^ u,v 
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Dropping the lower-order term that is linear in Yu^v now suggests the following approximation: 

Note that this expression is in a form that should provide self-correction in our estimate for Y^^v ■ If 
Yu^v is large compared with its expected trajectory then this term will be large and (as this term is 
negated in the expected one-step change) there is a 'drift' that pulls Y^^^ back toward the expected 
trajectory. 

Note further that our expression for 4:Pu^v + Ru,v + Su,v involves Yu and Y^. So, for this approach 
to be useful we will need some control on vertex degrees. We achieve this control using the same 
mechanism we are using to control co-degrees. Namely, we consider the expression in the numerator 
of the expected change in Yu, 

1 V- 



E[AYu\Ti] = -- y,Yu,, 



where 



{xy G E{i) :x,y£ N^} 



and track this quantity for every vertex u. This quantity should be correlated with Y^, that is, we 
should have 



Tu^p 
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We arrive at the ensemble of variables Q,Yu,v, Ru,v, Su,v, Pu,v,yu,Tu and Yu^v,w (it turns out that 
some control on the common neighborhoods of triples is needed for our martingale arguments), 
as illustrated in Figured! We are ready to state our main result, which establishes concentration 
for the entire ensemble of variables and in particular provides us with the level control over Q 
required to deduce Theorem [TJ Here and in what follows the notation A = B ziz C is short for 
Ae[B-C,B + C]. 



5 



Theorem 2.1. Letp{i,n) = 1 — 6i/n and set 

7 = 2-\/2, 7 = 7(n)=7-- and ^> = $(T, n) = e'^(^"P) log^ n . 

logn 

Then there exist absolute constants a, /3, 5, k, ^u, T > such that with high probability 

Q = nV/6 ± a^nV"^"^^^ (2-5) 

y„ = n/ ± an^/V^ (2-6) 

Ru,v = pYuAYu - Yu,.) ± /5n3/V+^$ (2.7) 

5n,. = pYuAYv - Yu,.) ± /3n3/V+^^ (2.8) 

Pu,,=pYlJ2 ± Sn^/^p^+''^ (2.9) 

Yu = np ± KU^'I'^p^^^^ (2.10) 

r„ = pYll2 ± (2.11) 

Yu,v,w = n/±2max|(np^log^n)^^^ , log^ n} (2.12) 

/or all u,v,w and throughout the time-steps < i < 1 — (6a^e^^ lo^_^^i/(4 27) ^2yg^ 

To deduce Theorem [U let p* be the lowest value of p for which Theorem 12.11 holds and note 
that 

/ 1 10 \ 1/(4-27) 

= ( sa^e^^i^^ ) = O fn-V(2v/2) i^g^/v^ (2.13) 



n 

since 4 — 27 = 2-v/2 + 0(1/ logn). In particular, p = p* satisfies n^p^/6 > a^n^p^^~^$^ since 
$ = e^(i~p) log^ n < log^ n. It thus follows that Q > w.h.p. due to (12. Sp and it remains to 
recover the number of edges corresponding to p*. Recalling (j2.3p we have 



\Eii)\ = (2) - 3z = Q - ^(1 - P'W = \ (nV - n) 

and the desired result follows from (|2.13p with room to spare in the power of the logarithmic factor. 

We prove Theorem 12.11 in the following section by applying martingale arguments to random 
variables that track the differences between the random variables we are interested in and the 
variables they should follow. Note that we establish some form of self-correction for every variable 
in this ensemble, with the notable exception of Y^- 

The authors suspect that the methods introduced in this paper can be further developed to 
achieve better high probability upper bounds on |£'(M)|. This might be achieved by expanding the 
ensemble of random variables (perhaps using ensembles of generalized extension-counting variables, 
which is the approach taken in the recent analysis of the H-hee process [6]). However, it seems 
that a nontrivial modification would be needed to prove the conjectured bound |£'(M)| = n^/'^+°W . 

For notational convenience we set 

A ' 



log^ n 



Note that while Thereom 12.11 applies, estimates (|2.6p - (j2.1ip and (j2.5p can each be written as a 
main term times (1 -|- 0(1/ log^ n)) = {1 + 0(A)). Throughout the paper we will use a convention 
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whereby all Greek letters are universal constants. We do not replace any of the constants (including 
the pivotal 7) with their actual values. This is done in the interest of understanding the role these 
constants play in the calculations; it turns out that these constants must be balanced in a fairly 
delicate way. We observe that these constants can take the actual values 

1 a X « 

a = - P = — 



2 '2 4-27 

(1 — 7)a a 

K = ^ u = - T = 20 . 

2 ^4 

The key conditions these constants must satisfy are ([3l9]l . ([322]), ([3:2i]) . ([3:28]l and ([5:29]) . The 
only condition on T is that it is sufficiently large. 



3 Proof of Theorem 12.1 

Define p* as in ([2T3]) and let i* = 1(1 — p*)n be the analogous round. Let Gi be the event that 
all estimates in Theorem 12.11 hold for the first i steps of the process. 

For each variable and each bound (upper and lower) in Theorem 12. II we define a critical interval. 
This interval has one end at the bound we are trying to maintain and the other end slightly closer to 
the expected trajectory of the random variable. If one of the estimates of Theorem 12.11 is violated 
then the corresponding random variable 'crosses' a critical interval. We bound the probability 
of each such event using a martingale argument, introducing a separate supermartingale for each 
variable and bound of interest and for each step in which the variable could enter the critical 
interval. Theorem 12 . 1 1 then follows from the union bound (note that the number of supermartingales 
we consider is bounded by a polynomial in n). We restrict our attention to these critical intervals 
because the expected one-step changes in our random variables each have a 'drift' term that pushes 
a wayward variable back toward the expected trajectory. By only considering the critical intervals 
we make full use of these terms: This is the mechanism we are using to establish self-correcting 
estimates. For an application of this idea in a setting with fewer variables, see [5]. 

Let the stopping time r be the minimum of i* and the smallest index i such that Qi does not 
hold. Consider an event £ of the form X{i) < x{t) for all i < i* where we assume that X{i) is a 
random variable and x{t) is not. Note that every bound in (I2.5p - (j2.12p can be written in this form; 
that is, the event {r = i*} can be written 

where |X| is polynomial in n and each event Eg is of the form X{i) < x{t) for all i < i*. For each such 
event £ we introduce a critical interval of the form l£{t) = {x{t) — w{t),x{t)) where w{t) = o{x{t)). 
Consider a fixed step iq, which we view as a step at which X{i) might enter the critical interval Ig. 
Set to = io/n^- Define the stopping time Tg^i^-^ to be the minimum of max{io,T} and the smallest 
i > io such that X{i) If. Note that if X(io) leito) then we have T£^ig = zq. Thus this stopping 
time is (formally) well-defined on the full probability space (n.b. we only make use of this stopping 
time when X(io) is in the critical interval and X{iQ — 1) is not). We now establish a bound B(i) 
on the one-step change in X{i) conditioned on Qi. This bound is far less than the width w{t) of 
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the critical interval. Given a particular event £ and starting step ig, we work with the sequence of 
random variables 



Note that if X{iQ — 1) is not in the critical interval (and < ''") then Z(io) < —w{tQ) + -B(io)- 
Therefore, in the event £^ there are steps io < j < i* such that Z{iQ) < —w{tQ) + -B(io) and 
> 0. However, our stopping time r stops all of these sequences as soon as any of our conditions 
(|2.5|) - (|2.12|) are violated. So, we have 



It remains to bound the probability of each event in this union. This is done for each of the bounds 
(j2.5p - (|2.12p in turn in Sections 13.11 - [3^71 In order to bound the probability of these events we will 
apply the following inequality due to Preedman [9] , which was originally stated for martingales yet 
its proof extends essentially unmodified to supermartingales. 

Theorem 3.1 ([9], Thm 1.6). Let (Sq, Si, . . .) be a supermartingale w.r.t. a filter (J-i). Suppose 
that — Si < B for all i, and write Vt = X]i=o ^ [('S'j+i — Si)'^ \ J^i] . Then for any s, -y > 



Our applications of this inequality will each have two parts: a careful calculation that establishes 
a martingale condition and a coarser argument that provides bounds on both the one-step changes 
and the second moment of the one-step changes of these variables. We emphasize that our carefully 
chosen stopping times allow us to assume that the event Qi holds throughout these calculations. 
This is henceforth assumed without further comment. 

3.1 Edges between a co-neighborhood and an exclusive neighborhood {Ru,v, Su,v) 

In this section we establish concentration for Ru.v,Su,v (by symmetry it suffices to consider Ru,v) 
throughout the process via a delicate analysis of the behavior near the tails of their distributions. 

We begin with a calculation of the one-step (conditional) expected change in Ru,v To simplify 
the exposition we will write xy G Ru,v to denote that the edge xy is counted by the variable Ru,v, 
and similarly for Su,v, Pu,v,Tu- There are 7 types of triangles that contribute to this expected 
change (see Figure [2] for an illustration): 

1. Triangles vxy where x S A''„,„ and y £ Ny\ N^^v and y ^ u. For a fixed x £ N^^v there are 
Yv,x — Yu,v,x — ^{uveE} such triangles and selection of one of these triangles moves x from Nu^v 
to Nu \ Nu^v and thereby decreases Ru,v by Yu^x — Yu,v,x — ^{uv<^E} (due to the neighbors of 
X in the exclusive neighborhood of u) and increases Ru,v by Y^^^^j. (due to the neighbors of x 




teX l<io<j* 



= U U {Z{n- Z{iG)>w{t^)-B{i^)}. 



£6X l<jo<«* 




8 



u 







(D 


1 




N 







i 



N,\Nu^,\{u] 



Fi gure 2: The types of triangles contributing to the change in Rn y. 

in Nu^v)- This results in a contribution of 

E Yy^x - yu,v,x - '^{uveE} ii _ _ \ _v \ 

Q \yJ^u,x -fit, II, X ^{uvGE}) J^u,v,x) 



xeNu,v 



E - 



xeNu,v 



^"1 ^"-"'^ (Y -2Y ) + o(- 



(3.1) 



where in the last equality we absorbed the indicator variables into an 0{p/n) term based on 
the estimates of Theorem 12.11 that hold up to the stopping time r. 

2. Triangles uxy where x S Nu^v and y £ Nu \ Nu^v, i-e. triangles given by xy G Ru,v The 
selection of such a triangle removes y from A^^j^ while moving x from A^^j^^^ to \ N^^^ . Noting 
that there are Yu^x — Yu^v,x — '^{uveE} edges of the form xy E Ru,v for each x £ Nu^v while 
there are Yu^v,y edges xy G Ru,v ^or each y G Nu \ Nu^v, selecting the triangle uxy decreases 
Ru,v by Yu^v,y + Yu^x — Yu^v,x — '^{uveE} ~ 1- Overall this results in a contribution of 



E 

x£Nu,v 



(Yu^x Yu^v,x) 

Q 



E 



Y^ /D 



Q 



n 



(3.2) 



yeNu\Nu,v 

where the indicator terms were again absorbed into the 0{p/n) term. 



3. Triangles vxy where xy S Pu,v Choosing such a triangle moves x,y from N^^v to A^^^ \ Nu^v 
and thus decreases Ru,v by {Yu,x - Yu,v,x - '^{uveE}) + (Y^^y - Yu,^,y - l{uv(iE}) and increases 
Ru,v by Yu^^^x + Yu,v,y - 2. This contributes 



^ ^ ^{Yu^X 2y^ 5; a; + Yu^y 25^ ,y y) + O | j 



xyePu.v 



xeNu,v 



Q 



(3.3) 
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4. Triangles uxy where xy G Pu,v Selection of such a triangle moves x, y from N^^^ to \ Nu^v 
and so decreases Ru,v by {Yu^r, - Yu,v,x - '^{uveE}) + {Yu,y - yu,v,y - l{™e£;}) (similar to the 
previous case only without introducing any new edges to Ru,v)- This amounts to 



^ \''-U,X ^u,v,x) "-^ \ 



xyePu,v 

= - E %^(^-,- - Yu,v,x) + O f ^) . (3.4) 

xeNu, 

5. Triangles uxy where xy is in the set of edges induced by Nu \ Nu^v , i-e. xy G Tu \ {Ru,v U Pu,v)- 
Each such triangle decreases Ru,v by Yu^v,x + This results in a contribution of 

xy6r„\(ii;„,„UP„,^) ^ x(^Nu\Nu,v 



-Hu.m E n = n = ^\^)- ^^-^^ 



6. Triangles uvx for x G A'lj^^ contribute 

7. Triangles xytj; where x S A^m,„ and ?/ G Nu \ N^^v and w ^ u. We note in passing that this 
is the only type of triangle whose selection impacts Ru,v while changing neither Y^ — Y^^v 
nor Yu^v For a fixed xy G Ru,v there are Y^^y — 1 such triangles and this results in a total 
contribution of 

to the expected change in Ru,v If we simply apply our bounds on Yx^y to the expression (13. 7p 
we would get an undesirable error. A more careful approach will give better bounds. Namely, 
we rewrite this expression by summing over x G N^^^ and using our error bounds on Su,x and 
Pu,x- This should give a better estimate as Su,x and Pu.x 'aggregate' many terms Y^^y and 
should therefore produce better error bounds, and indeed this proves to be a crucial choice. 

There are two types of triangle xyw in this category depending on whether or not uw is an 
edge. Consider first the triangles where uw is not an edge. Observe that an edge yw which is 
in Su,x due to having y G Nu.x and w G Nx \ Nu,x determines a triangle of this type so long 
as y V and y ^ N^. The requirement y ^ v excluded Su,v edges from being yw candidates 
provided that uv G E, whereas edges excluded due to y £ correspond to having x,y £ N^^v 
with xy £ E and w G Nx^y \ {Nu^x,y U {u}). Therefore, the number of triangles of this type is 

^ ^ Su,x — ^{uv&E}Su,v — 2 \Yx^y — 2 — Yu^x,y + l{tn;e-E}] 

xeNu,v xyePu,v 

= E [pyu,x{yx-Yu,x)±/3n^^V+"'<^]- E ^[yx,y-Yu,x,y] + 0{nY) 

xeNu,v xyePu^v 

= 5^ pYu,x{Yx-Yu,x)- E 2[y,,j,-y,,,J±(l + 0(A))/3n5/y+^$, (3.8) 

xeNu,v xyePu.v 
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where in the last equality we absorbed the O(n^p^) error term using the fact that A = log ^ n 
is 0(n^/^p^) (with room to spare). 

We now turn to treat triangles where uw G that is, triangles xyw which either have 
w,y € Nu \ Nu^v and x G N^^v or alternatively y G Nu \ N^^v and ui, x G N^^v Note that the 
removal of such a triangle will decrease Ru^v by 2. The effect of removing these triangles can 
be counted by ^^^.g^y^ ^ Pu,x using a delicate double-count: In the case where w G iV„ \ Nu,v 
the aforementioned sum counts the triangle xyw once, whereas this triangle is counted twice 
if u; G Nu^v Furthermore, this sum includes pairs yw G Pu^x with w,y £ Nu^v, which do not 
affect Ru^v at all. Weighing in the number of triangles xyw as above together with their effect 
on Ru^v now totals to taking 2Pu^x, negating 1 per triangle with w G N^^jj and another 1 per 
triangle with w G Nu,v,x (the last case is counted twice by the edge yw with w,y £ N^^v and 
hence cancels away the contribution of this triangle) . Altogether this sums to 

2{Pu,x - l{uvmyu,v,x) - Yl {Wu,x,w\iNu,v^{v})\+\Nu,v,x,n,\) 

XGNu,v 



E 



'2{Pu,x ~ '^{uv^Ey^u,v,x') ~ 



E 



(^^,2;,«j 



'-{uvi^E 



}) 



X&Nu,v 



Y ± 25n3/V+^$ - Y ^™ + 0(V 



weNu,v,: 



E 

xeNu,v 



weNu,v,x 

Putting ()3.8p and ()3.9p together, taking advantage of the fact that trivially 

Y^ 2yn,x,y = Yl Yl ^u^x^w ) 
xyePu,v xGNu,v weNu,v,x 

we see that the expression in (j3.7p is equal to —l/Q times 

Y {pY,,,,Y,±36n'/V^^^^-2(^ Y ± (1 + 0(A))/3n5/y+^$ 

xeNu,v xyePu,v 

= Y [p^uA^^P ± «^'/V"'$) - Yu,,,x{np^ ± an^/V^-) 

x£Nu,v 

± [(1 + 0(A))/3n^/V+^^ + 45n5/^/+^$ 



(3.9) 



Ru,vnp ± 



[1 + 0(A))Kn^/ V"^"^^ + 6an^/V^^^ max{np^ log^ n} 
+(1 + 0(A))/3n5/V+^^ + 45n^/V+^$ 



(3.10) 



where we used the fact that R^^v = XlxeA^u „(^">a;~^f>'^>2:"'"^{M«^e-E}) ^'^d absorbed the indicator 
variables (whose total contribution is at most np'^Y^^v = 0{n?p'^)) into our final error term, as 
well as the inequality Yu^v,x ^ 3max{np'^, log^ n} which follows from (12.120 by monotonicity. 

Now that we have an expression (albeit in 7 parts) for the expected change in Ru,v7 we are ready 
to consider the expected change in Ru^v relative to its expected trajectory. Define 
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and consider E[AX | Write 

A [p{Yu - = p(i + 1) • A [(y„ - F,,,)y„,,] + Ap • (y„(i) - y„,.(i))^n,.(0 • (3.ii) 

We will see that the expected change in Ru,v that comes from triangles of types 1-6 will balance 
off with the first term in ()3.1ip , while the expected change in Ru^v due to triangles of type 7 will 
be balanced with the second term. 

Collecting (j3.ip - ()3.6p . the total contribution to E[Ai?M^t, | Fi\ from triangles of types 1-6 equals 

^ ^ [^,a;-^,a; ~l" ^,a; ~ Yu^xYu,v,x '^Yy^xYu^Vjx'] ^ ^ ^ Yy_^xYu,v,x ~l~ O ^ ^ . (3.12) 

X(iNu,v x£Nn\Nu,v 

Furthermore, we can analyze the change in H = Yu,v(Xu — Yu,v) by considering the following 3 cases: 

(i) Selecting a triangle uxy for x,y & Nu \ Nu,v- the co-neighborhood of u,v is unchanged while 
(Yu — Yu^v) is decremented by 2, altogether giving AH = —2Yu^v- 

(ii) Selecting a triangle vxy for x G N^^v aud y ^ u: the co-neighborhood loses x and in addition 
loses y if y G iV^u^t,, while 1^ remains unchanged. The change in H is thus 2Yu^v — Y^ — 1 in 
the former case and 41^^^, — 21^ — 4 in the latter one. 

(iii) Selecting a triangle uxy for x G N^^v- If y £ the co-neighborhood loses 2 vertices while 
A(y„ - Yu,v) = 0, yielding AH = -2(y„ - y„,„). Similarly, if y = v then AH = -(y„ - y„,„). 
Otherwise, y G \ Nu^v and both Yu^v and (y„ — decrease by 1, hence AH = —Y^ + 1. 

Altogether, we can obtain the factors of 2 in Item ^ and in the case y G Nu,v of Items . (jui]) 
automatically by symmetry when summing over x as follows: 

E[A[y„,, (y„ - y„,,)] | f,] = -y.,, -^"-^ " ^"•''•^ 



(y.-2y„,.-o(i)) ^"'^ ^^"""''^ 



All the triple-degree terms cancel out and we can sum all the 0(l)-terms and rewrite the above as 

^ ^ [y!i,3;(yu Y^^y^ + Yy x{Yu — ^ ^ Yy xYu,v ~^ ^ ( ) ■ 

Notice that when multiplying the above by p the error term becomes an additive 0{p/n) while the 
main terms are of order 0{p^). As such, the same estimate holds for the result of multiplying the 
above by p{i + 1) (which differs from p{i) by an additive 0(n~^) error and thus introduces an extra 
0{{p/n)'^) error term). We can now combine this with the change in Ru^v given in (j3.12p to get 
that the contribution to E[AX | Ti] from triangles of types 1-6 and the first term in (j3.1ip is 

^ ^ [ {Yu.xYv^x ~\~ Yy^ ,j. Yu xYu^v,x '^Yy xYu^v,x^ ~^ pYu^x(Yu Yn -n^ ~\~ pYy xiYii 2yj^ ,y)j 
xeNu.v 
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xeNu\Nu,v 

which is equal to 



p 



+ 

x€Nu\Nu.v xeNu,v 



(3.13) 



The key observation here is that the first and second summations in (|3.13|) feature the random 
variable X itself, a fact which our self-correction argument for X hinges on. Namely, by definition 
of Ru,v we have Z]xeAf„,„ [^u,x - yu,v,x - '^{uveE} - p(Xu - = Ru,v - pyu,v{Yu - Yu^y) = X 

and similarly "^^eNuXNu v^^'^'^'^ ~ P-^",^] ~ -^^ hence we can use our (1 + 0(A))-approximation for 
the factors {Yu^x + in the summations in (j3.13|) and obtain that the contribution to 

E[AX I J^i] from triangles of types 1-6 and the first term in (j3.1ip is 



12 + 0(A)^ _ 6 + 0(A) ^ 6 + 0(A) 



18 + 0(A) ^^ ^ 12 + 0(A) ^ 



n?p 



n^p 



E 

xi^Nu.v 



+ o(^ 

n 



(3.14) 



We now turn to the triangles of type 7. As we noted above, we balance the term (13.10p with the 
second term in (j3.1ip . i.e. the expected change in pY^^viYu — Y^^v) due to the change in p (which 
deterministically decreases by 6/re^). The total contribution to E[AX | Fi] from these terms is 



Ru,vnp^ 6 



Q 



~\~ <)Yu v{Yu Yy^ lb 



(1 + 0(A)) (/? + K)n^/V^^^ + 45n5/V+^$ 
-|-6an^/^p^'^'^<I> maxjnp^, log^ n} 



6-bO(A) r »^+^$ 7)2+7$ 



n^p 



n 



n 



1/2 



(3.15) 



where we used the fact that the term corresponding to taking log^ n in maxjnp^, log^ n} has order 
72~3p7-i ^Qg5 = 0(log^n • n~^/'^p^~^^(^) as long as p > n~^/2log^n (which easily holds for large 
enough n) and hence is negligible compared to the 0(A)-error in the term involving 13 + k. 
Putting everything together, the combination of (j3.14p . (j3.15p gives 



E [AX I Ti] 



24-bO(A) r n^+^^ 

^-X±(6 + 0(A)) {i3 + K)^-jj^ + i6a + 65f 



T1?P 



n 



n 



1/2 



(3.16) 



and we are ready to establish the concentration of Ru,v via a martingale argument. As outlined 
above, we introduce two critical intervals for the random variable X, corresponding to the upper 
bound and lower bound on Ru,v These intervals have one endpoint at the bound we are trying to 
establish and the other somewhat closer to zero (corresponding to the expected trajectory of X). 
For the variable Ru,v to violate Eq. (|2.7p it must be that X crosses one of the critical intervals. 
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Our critical interval for the upper bound on Ru,v is 

= (/Jn^/ V"^"^^ , /3n3/V^^^) where /? = (l - log"^ n) j3 . (3.17) 



Suppose that X{iQ) enters Ir for the first time at round io (within the time range covered by 
Theorem 12. ip and define the stopping time tr = min{i > : X{i) < ^n^/^p^+^^>}, i.e. the first 
time past that point at which X exits the interval through its lower endpoint. We claim that 
Z{i A tr) is a supermartingale, where 

Z{i) = X{i) - /3n^/'^p'^+'^^ for i > io ■ 

To see this, write t = i/n? according to which p = 1 — 6t and <I> = e^^*log^?i, and note that for 
any 7 > the second derivative of f{t) = e^^*(l — 6^)^'^'^ is uniformly bounded in [0, |], hence 

f[t + = f{t) + (6Te6^*(l - 6^)2+^' - 6(2 + j)e^^\l - 6t)^+"') n'^ + 0(?i-^) . (3.18) 



This provides an estimate for A[/3n'^/^p^+^<I>] between Z{i + 1) and Z{i). At the same time, for any 
X{i) satisfying (|3.17p we can plug in the lower bound this gives for X in (|3.16p and obtain 



E[AZ \J^i,TR>i]< -{24 + 0(A)) /3^— ^ + (6 + 0(A)) (/3 + 



1+7 <|) „...x,„ ,p^+^<I) 

jl-L/^ ' ' ' ' ' ' 77,1/2 



K 



+ (36 + 0(A)) (a + S)^^ + (1 + 0(n"2)) 



r)H-7$ 752+7$ 

6(2 + 7)/3%^-6T/3^-^ 



4/3 + 3/3 + K + 7/3 + 0(A)) ^ + ( - + 6a + 65 + 0(A)) -^-^ 



By the definition of /3, 7 we have —4/3 + 3/3 + k + 7/3 = — /3 + k + 7/3 — 2/3/ log n as the difference 
between 7 and 7 negates the difference between /3 and /3 as well as the 0(A) term, hence we can 
conclude that if n is large enough then Z{i A tr) is a supermartingale provided that 

(l-7)/3>« (3.19) 

and T is sufficiently large, namely T > 6(a + (5)//3. 

To apply Freedman's inequality we need to obtain bounds on \ AZ\ and E[(AZ)^ | J^i]. Recall 
that AZ = ARu,v - A[pH] - A[f3n^/^p'^+''<^] where H = (y„ - Yu,v)Yu,v By (IXTTD and the fact 
that Ap = we have A\pE] = pAE + 0{p^), while A[/3n3/ V"^"^^] = 0{n-^/'^p^-^^^) = o(l). 

Hence, AZ will be dominated by AR^^v — ^b^]- There are four cases to consider here: 

(i) Choosing a triangle that includes u or u and some vertex x G Nu^v (triangles of types 1-4,6 in 
the study of K[ARu^v \ J'i] above): There are 0{v?'p^) such triangles and selecting one of them 
affects both Ru,v and H. As we next specify, the principle terms in these changes are identical 
and so \AZ\ is bounded by the error terms in our approximations. Indeed, going back to the 
analysis of the triangle types as well that of AH we recall the various triangle types satisfied: 



• Type 1 

• Type 2 

• Type 3 

• Type 4 



ARu,v = -Yu,^ + 2Yu,^,^ + 0(1) vs. AH = -Y^ + 2y„,^ + 0(1). 

ARu,v = -Yu,x + {Yu,v,x - Yu,v,y) + 0(1) vs. AH = -Yu + 0(1). 

ARu,, = 2{Yu,,,x + Yu,,,y) - {Yu,x + Yu,y) + 0(1) vs. AH = AYu,, - 2Y^ + 0(1). 

ARujV — (Yu,v,x ~l~ ^^,i),j/) (Yn x ~l~ Yu y^ + 0(1) ^s. A,^ — 2Yu^v 21^. 
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• Type 6: ARu,v = -Ru,v vs. AH = Yu^v - Yu- 
In all of the above cases the main terms of Ai?u^„ cancel those of pAH at the cost of an 
0(ni/V iog2 

n)-eiTor due to the approximations of Y^^v and Yu (dominating all other errors). 

(ii) Choosing a triangle that includes u but no vertex in N^^v (triangles of type 5 above): there are 
O(n^p^) such triangles and each corresponds to ARu^v = — {Yu,v,x + 'Yu,v,y) vs. AH = —2Yu^v 
Hence, in this case \AZ\ = 0{\/ np^ log^ n + log^n), with the dominant error-term due to 
triple co-degrees (exceeding that in pYu^v since 7 > |). 

(iii) Choosing triangles that include an edge of Ru,v but contain neither u nor v (type 7 triangles): 
Each of these 0{n^p^) triangles corresponds to \AZ\ = 0(1) as AR^^v = —1 and AH = 0. 

(iv) Choosing any other triangle: as R^^v and H are both unchanged, these triangles can modify 
Z by at most 0{p^) due to the additive 6/n^ change in p. 

The L°° bound on AZ is clearly dominated by rounds of the first sort and \ AZ\ = 0{n^^^p^ log^ n). 
For an bound notice that the 4 round types contribute 0{p^~^'^^ log^ n), 0(p^log^n + " ), 
0{p^) and O(p^) respectively to E[{AZ)'^ \ J^]. As 7 < 1 we have log^ n = 0(^^+2"^ log^ n) 
while the fact that p > p* (which has order n~^^^^~'^^^^°^^^ as was seen in Eq. (j2.13p ) implies that 
= 0(^^+2^-) since np^+'^'f > n(3-47)/(4-27)-o(i) > ^1/5 ^j, ig^j,gg enough n. Altogether it follows 
that E[(AZ)2 I Ti] = 0(y+27iog4^). Clearly the L°° and the bounds on AZ also hold in the 
conditional space given tr > i. 

Recall that we are interested in Z{i) starting at time iq, i.e. immediately after X enters the 
critical interval Jr. Let pq = p{io) = 1 — Gio/n^ and observe that our bound on \AZ\ guarantees 
that < X{io) - /3n3/2(po)27$ < 0{n^/'^{poV log^ n). Hence, 

Z{io) < 0-f^)7i'/Hpof-^^<^ + O{7^'/\po)nog'7^) < -i/3n3/2(po)2+nogn, 

where the final factor of | readily cancels the 0{n^^^ (pq)^ log^ n)-term for sufficiently large n since 
n(po)^log"^n > ni-°(i)(p*)2 > n^/^. We now apply Freedman's inequality (Theorem 13. ip to the 
supermartingale Sj = Z(^{iQ + j) A r/j) while noting that the above analysis implies that 

So<-i/3n3/2(po)2+nogn, 
max \ Sj-^i — Sj\ = O (^n^^'^{poy log^ nj , 

J;e[(5,+i - S,r I J-j] =o(log'nY,{p{i))'+'A = O [n'ipor'nog' n) , 

j ^ i>io 

where J-^- = We deduce that for some fixed constant c > 0, 

P(U,{5, > 0}) < exp|-c f^^f^'^^^j'Z ^1 =exp{-m(po)'log-'n} , 

VJlJ- n2(po)2+27log5n + n2(po)2+27log3nj ^ ^' 

which is sufficiently small to afford a union bound over all u,v and time steps iq, hence w.h.p. X 
never crosses the critical interval Ir and in particular X{i) < /3n^/2p2+7<|) ^^j. ^ ^^^^ Exactly 
the same argument shows that w.h.p. X{i) > — /3n'^/2p2+7<j) for all u,v,i by considering the lower 
interval ( - /3n^/2p2+7<j) ^ -/3n^/2p2+7<^) and analyzing the variable Z{i) = -X{i) - /3n^/2p2+7<|). 
This completes the proof of Eq. l^l\\.dM. 
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3.2 Edges within the co-neighborhood {Pu,v) 

Our analysis of the variable Pu,v (as well as that of the remaining variables throughout this section) 
will follow the exact same arguments used in Section 13.11 to control Ru,v , namely we will study the 
variable X = P^^v — P^uvi establish L^,L?',L°° bounds for its change per round and deduce the 
result from an appropriate martingale argument. This being the case, we will repeat the arguments 
only briefly while keeping the main focus on deriving the required estimates on A.X. 

Observe that choosing any triangle incident to some x G N^^v and either u or u (or both) removes 
X from the co-neighborhood and thus eliminates |A^M,t,,a:| edges from the tally of Pu,v Removing 
any other triangle xyz (i.e. one which includes neither u nor v) that includes 1 < s < 3 edges from 
the induced subgraph on N^^^ will decrease the count in P„ „ by these s edges. Overall we get that 



Y — 2 



= - E ^^^^^^ w - E % + o (t^) ' (3-20) 

where the term —2Pu^v/Q is due to the fact that any xy £ E with x,y £ N^^^ (a total of Pu,v 
edges) is double-counted in the first summation upon selecting one of the triangles uxy or vxy. 
Next define 

X{i)=PuAi)-pYlJ2, 

and note that since Yu^v{i) — 2 < Yu^i,{i -|- 1) < yu,v{i) we can write A(Y^^) as {2Yu^v + 0{l))AYu^v 
and using the fact that p{i -|- 1) = p{i) — 6n^^ as well as Eq. (12. ip we deduce that 

E[A(py,y2)|^,]=--^^ + i(p--^)(2y., + o(i))(- E "^""^"^^''"^'"0 

= -^-pY^,, E • (3.21) 

Combining (|3.2U|) . (|3.21|) and noticing that by definition "^xeNu v^^'^''"'^ — pYu,v) = 2X now gives 



xeNu.v xyePu,v 



12 + 0(A) ^ _ V±anVV£ ^Y^, ( 



^ X(iNn,v ^ ^ 

^ 30 + 0(A) ^ (3 + Q(A))a/+^^ 
n^p ' 

where we used the approximation that all co-degrees, Q and Pu^v are (l + 0(A))np^, {^ + 0{K))n?p^ 
and (^ + 0{K))'n?p^ respectively at time i by our hypothesis. 

Define 5 = {\ — log~^ n)5 and let assume that in some round io the variable X(io) enters the 
upper critical interval given by 

Ip = ('5n^/V+T'$ , 5n^/V+^$) where 5 = {l - log"^ n) 6 . 
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Let Tp denote the stopping time min{z > io : X(i) < Sn^^'^p'^^^^} and consider the random variable 

Z{i) = X{i) - (5n3/V^^^ for i > io . 

Our upper bound on E [AX \ and the fact |A[(^n3/2p3+7$]| < (6 + 0(n-2))(3 + 7)<5n-i/2p2+7$ 
(derived by the exact same argument as in Eq. (j3.18p in Section I3.ip now imply that 



E[AZ \J'i,Tp>i]<3 -10(5 + a + (6 + 2j)6 + 0(A) 



n 



1/2 



Since the term 2(7 — 7) = 12/ log n cancels 10(5 — 6) = 10/ log n as well as the 0(log ^ n) term 
represented by A, it follows that Z{iQ + j A Tp) is a supermartingale provided that 

(4-27)5 > a. (3.22) 

To obtain L°° and bounds for the one-step changes in Z argue as follows: 

(i) If we choose a triangle that includes an edge joining {u, v} with N^^v this then modifies 
both Pu^v and pY^^/2 yet by the same leading order terms and hence changes Z by at most 
O ( \/ np^ log^ n + log^ n). The number of these triangles is 0{'n?p'^) as the dominant term here 
corresponds to edges in Ru,v, Su,v Thus, the probability of such an event is at most 0{p/n). 

(ii) Selecting a triangle that includes an edge in Nu^v but does not include u 01 v (and hence does 
not change Yu,v) modifies Z by 0(1). The asymptotics of Pu,v and the co-degrees imply that 
there are 0{n^p^) such triangles and hence the probability of such an event is O(p^). 

(iii) All other triangles do not change Pu,v,Yu^v and hence change Z (via the change in p) by 0{p'^). 
Overahweget \AZ\ = 0(nVV/^ log^^^ n+log^ n) andE[(AZ)2 | J^,] = 0(p^ log^ n+(p/n) log^° n) . 
Letting pQ = p{io) = 1 — Qio/n'^, the definition of io and our bound on |AZ| guarantee that 
< X{io) - (5713/2(^0)3+7$ < O{n^/^{po)'flog^n). Hence, 

Z{io) < {5 - 6)n^/\pof'^"'^ + 0{n^/\pof/Hog^/^ n + log^ n) < -l5n^/'^{pof+^ logn , 

where the final factor of ^ cancels the O(^n^^^{po)^^'^ log^^^ n + log^ n)-term for large enough n since 

„p3/2+7 log-3/2 n > „(f-37)/(4-27)-o(l) > ^1/4 n^/^p^+nog-^ n > n(3-47)/(4-27)-o(l) > ^1/4 

for all p > p*. We conclude that 5^ = Z[{iQ + j) A rp) is a supermartingale with 

So<-^5n^/\pof+nogn, 
max \Sj+i - Sj\ =0 n{pQ)^ log^ n + log^ , 

5^E[(5j+i-5j)2 I J^.] =0{n\pof log' n + nipof log'' n) , 
j 

where J-j = J'ic+j- Recalling that 7 > | and that in addition n~^/'^logn = o{p'^~^) for all p > p* 
since {p*)'^~'^ > cn~^^'^log' n by (|2.13|) . we see that 

|5o| maxISj+i - Sj\ = O (^ii^ (po)^^^ log^ n + n^pof+'^ log^ = o {n"^ (pof log' n) 

and therefore Theorem 13.11 gives that for some fixed c > 0, 

p,,ro>nl^^ / n3(po)6+27iog2n ) [ n2(po)^+27 

F {Uj {Sj > 0}) < exp <^ -c g ^0 ^ = exp " 

L n^(po) log n + n(po) log n) 



n{po)^ log n + log n , 
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Observe that if po > n-V3 then the last expression is at most 

exp{-ni-°W(po)'+'^} < exp{-ni(i-^)-°W} = exp {-ni(v^-i)-°W} < e""'^' 
for large enough n, and otherwise (recalling that po > > it is most 

exp {-n2-°(i)(po)'+'^} < exp {-n^^-^-^^i)} = exp j-n^'^v^-^d) } < e""'"' . 

In either case, the probability of the event Uj{5j > 0} is far smaller than any polynomial in n, hence 
a union bound establishes that w.h.p. X{i) < Sn^^^p^'^^^ for all u,v and i. A similar argument 
shows that w.h.p. X{i) > —Sn^^'^p^'^^^ for all u,v and i. This completes the proof of Eq. (j2.9p . 



3.3 Co-degrees (Fui,) 

Following the arguments in the previous sections we will establish L^, L^, L°° bounds on the variable 



X = Y„ 



np 



As p decreases by with each step, Eq. (j2.ip and (j2.4p show that 



E[AX I JFj] 



E 



+ - (2p - 6/ 



n 



n 



Q n 



n^p 



where the last error term absorbed the indicators and 0{n ^) from the first expression. Substituting 
our estimates (|2.7p . (|2.8p . (|2.9p for Ru,v, Su,vi Pu,v implies that this is equal to 

V / n yw^p J 

and using the estimate (j2.10p for Yu and that Q = (1 + 0(A))n'^p'^ we can conclude that 



n 



3/2 



n 



3/2 



+ 0( — 

n'^p 



Crucially, we did not approximate Q in the first expression with a (1 + 0(A)) correction factor as 
this would incur an error that would be too large to handle. Instead, there we apply Eq. ()2.5p and 
the fact that Yu^v = X + np^ to obtain that 



p(2np±2Kni/V^i$)y„,^ 12 + 0(A) 



Q 



n^p 



X + 



2v?p'^ 



1 j^3p3 _|_ Q,2j^2p27 1^2 



T^±(12 + 0(A)) 



Kp 



12 + 0(A)^^12g/ 

n \ 



rfip 



p27-4^2 



± (12 + 0(A)) 



j^3/2 



n J J ■ ■ re^/^ 

Combining this with the above estimate for E[AX | J^j], the term Vlpju vanishes and we get that 

'p27-3^2\ 



E[AX I T.^ = + ± (12 + 0(A)) 



n^p 



2bp'^^ ^ iip^-^^ ^ Kp"^-^^ 



n 



3/2 



n 



3/2 



n 



3/2 



+ 
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where the 0(l/(n^p)) error term was absorbed mto the O (n ^p^'^' -term since 7 < 1 and so 
l/(n^p) = o{p'^"'~^^'^ /v?). Furthermore, we claim that one may now omit this latter error-term 
altogether as it is negligible compared to the error-term of 0(A) in the terms involving /?, k. Indeed, 
keeping in mind that <I> and are each of order log^ n, we have 

„2 ^3/2 ^1/2^2-7-^^'^^ ^^(^.)4-27 n^/^ log^ n 

(3.23) 

Assume now that X{i) enters the upper critical interval defined by 



Iy = ian^^'^p'^^ , an^^'^p"'^] where a = (l — log ^ nj 



a . 



That is, suppose that is the first round at which X(i) > an^^'^p^^ and define the stopping time 
Ty = min{i > iq : X{i) < an^^'^p^^}. As before, consider 

Z{i} = X{i) - an^/'^p''^ . 

By the same argument of (j3.18p we have 



A 



an^^'^p^^ 



6 + 0(n-2)) (Tp"' - 7p^"M an-^/2$ 



and combined this with the above upper bound on E [AX \ J^i] establishes that 

E[AZ\Ti,TY>i]<[-2a + 2^ + 2K + ja + 0(A)] ^^3^^ + [4,5 - Ta + 0(A)] 

<[-(2-7)a + 2/3 + 2K]^^+[45-Ta + 0(A)]-|^, 

where we used the fact that 7 — 7 = 6/ log n to absorb both 2(a — d) = 2/ log n and the 0(A)-term 
for large n. In particular, Sj = Z{{iQ + j) A ry) is indeed a supermartingale so long as 

(2 - 7)a > 2/3 + 2k (3.24) 

and T sufficiently large, namely T > 45 /a. 

It remains to bound AZ in and L^. Here there are 2 types of rounds: ones in which we 
choose a triangle that involves u or t; and a vertex in Y^^v (there are 0{n'^p'^) such triangles) and 
ones where we choose any other triangle, in which case Y^^v is unchanged. The former event has 
probability 0{p/n) and leads to an 0(1) change in Z while the latter gives a variation in Z of order 
0{p/n) due to the — 6/?i^ change in p. Therefore, \AZ\ = 0(1) and E[(AZ)^ | Fi] = 0{p/n). 

Let pq = p{io) = 1 — Qio/n'^. By the definition of iq and the fact that \ AZ\ = 0(1), 

Z{io) < {a-a)n^/^{poy^ + 0{l) < -^an^^^ipoy logn 

(the last inequality holds for large n as the final expression clearly tends to 00 with n), and therefore 
the supermartingale Sj = Z[{io + j) A rp) satisfies 

So < -iani/2(po)^logn, 
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max — Sj\ = 0(1) , 
j 

j 

where = J^i^+j- Since n^/^(po)^~^ > clog^n for po ^ due to Eq. (|2.13p it follows that 
I Sol maxj l-Sj+i — -Sjl = O(n(po)^) and so Theorem 13.11 yields that for some fixed c > 

F(U,{5, > 0}) < exp{-c!^(^^2g^} = {-c(po)^-^log^n} < e"--^", 

which is sufficiently small to show that w.h.p. X{i) < an^^'^p'^^ for all u,v and i. The same 
argument handles the analogous symmetric case of the critical interval 

and shows that w.h.p. X{i) > —an^^'^p^^ for all u,v and i. This concludes the proof of Eq. (12. 6p . 



3.4 Edges within a neighborhood (T^) 

The number of edges in the subgraph induced by the neighborhood of u can change in two ways: 
Either a vertex is removed from (due to selecting a triangle of the form uxy with x £ N^) thereby 
decrementing Tu by all edges incident to it in this induced subgraph, or Nu remains unchanged 
(upon selecting a triangle that does not include u) and yet some of its inner edges are removed. 

The former case will be handled by directly summing over x E A^^^ , noting there are Yux triangles 
of the form uxy whereas the vertex x is incident to Y^x edges counted in Tu- The latter case requires 
a more delicate treatment, similar to the one used to study Ru,v in Section [3.11 Indeed, the naive 
approach would be to sum over edges counted by Tu, i.e. xy £ E with x,y £ Nu, as each of 
these would decrease T„ by 1 upon selecting one of the Y^^y triangles incident to it. However, the 
cumulative error in this approach (summing the co-degree errors for each edge in T„) would be 
quite substantial as it completely ignores the effect of averaging the co-degrees over T^. To take 
advantage of this point we will use our estimates for the random variables Su,x and Pu,xj which 
incorporate this averaging effect. Namely, we have 

E[ATJ J-,] = -l ^ y„2^ + ^-i ^ Pu,x-:^ Yl (3.25) 

where the first two terms account for triangles of the form uxy (the correction term Tu/Q is due 
to the double-count hy x,y £ Nu), the third term accounts for triangles of the form xyz where 
x,y,z £ Nu and the final term accounts for triangles of the form xyz where x,y £ Nu and z Nu- 

We evaluate the last two terms in (j3.25p using the bounds (|2.8p . (|2.9p for Su,v and Pu,v to get 
~ i2Pu,x + Su,x) = -^ Yl (p>;%+P>"x,«(^x-^x,.)±/3n3/V+^$±2<5n3/V"'^<J>) 

Y {yxYx,u) ± (1 + 0(A)) 



^ p_ 

2Q 

^ x&Nu 

We now need the following simple lemma. 



(3.26) 
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Lemma 3.2. Let ai, . . . , Um G IK and suppose that \a — ai\ < e for all i and some a G M. Then 



m 



+ 4me^ 



m 



Proof. The lower bound is due to Cauchy-Schwarz. For the upper bound fix 6 = ^ Oj and 
note that the assumption on the Oj's implies that |6 — ajj < 2e for all i. Observe that the convex 
function J2i of achieves its maximum over the convex set = bm} H (p|^ {\b — ai\ < 2e} ) at 

an extremal point where Oj = 6 it 2e for all i. Hence, Oj = 6 — 2e for [m/2j indices, Oi = b + 2e for 
[m/2\ indices and if m is odd there is a single = b. In particular 



< [m/2\{b - 2ef + [m/2\{b + 2e)^ + b'^l{m = i{mod2)} < mb"^ + 4m£^ 



□ 



i=l 



Revisiting Eq. ()3.25p . apply the above lemma with m = 1^ = (1 + 0{A))np and e = an^/^p'^^ 
to estimate the first sum there (noting that Xl^-g^^ Yux = 2Tu) and then use (I3.26P to get 



E [ATu I J^i] 



- 



Q\Yu 
±(1 + 0(A)) 



n 



1/2 



n 



1/2 



+ 0(l/n) 



(3.27) 



As usual set 

X{i) =Tu- pYll2 

and consider A(pl^/2). Observe that changes if and only if the triangle selected is of the 
form uxy with x,y £ N^, in which case it decreases by 2. Hence, E[A1^ | J^i] = —2Tu/Q and 
A(yjf ) = {2Yu — 2)A1^, and putting these together we get 

Combining this estimate with (I3.27P and the bound (I2.10|) for Y^ gives 



E[AX I J^i] 



4r, 



QYu 
±(1 + 0(A)) 



y2 



2q[ ' 



■ i^ AT / 



3y;2 



n 



1/2 



n 



1/2 



n 



+ 0(^ 

n 



By Eq. (|3.23p we have n ^p"^^ = o(An ^/^p'^^) for all p > p*, thus the above expression 
involving can be absorbed into the 0(A) error-term of the expression involving /3. Furthermore, 
since 4T„/(Qy„) = (12 + 0{A))/{n'^p) and 



2Q V 



we can conclude that 



n 



n'^p 



n 



/2 



E[AX I J-,] 



i^±^X±(l + 0(A)) 



3(/3 + k)p'^<^ ^ 6V+^$ 



n 



1/2 



n 



1/2 
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Now we consider the upper critical interval for given by 

Irp = ( , /in^/^p^+'^<^ ) where fi = (l - log~^ n) /i . 



and as before let iq be the first round in which X(i) > fin^^'^p^'^^^, define the stopping time 
TT = min{i > io : X{i) < /tn^/^p^"'"^$} and consider 

Z{i) = X{i) - fin^/Y+"'<l> . 

Exactly the argument of (I3.18P gives 



A 



3/y +7^1 = (6 + 0(n-2)) (Tpi+^ - (1 + 7)p^) nn^^/^<^ 



and together with the above upper bound on E [AX \ we get 

E[AZ\J^i,TT>i]<[-6fi + (3 + K + i2 + 27)/i + 0(A)] +[6-Tfi + 0(A)] -^^^ 



<[-(4-27)M + /3 + «]5^+[5-T^ + 0(A)]^^, 

where we used the fact that 7 — 7 = 6/ log n to absorb both 4{fi — fi) =4/ log n and the 0(A)-term 
for large n. In particular, 5*^ = Z{{iQ + j) A ry) is indeed a supermartingale provided that 

(4 - 27)// >/3 + K (3.28) 

and T is sufficiently large, namely T > 6/fj,. 

Having established an bound on AS it remains to consider the corresponding L^, bounds. 
If we choose a triangle of the form uxy, an event that has probability 0(l/n), then T„ decrease 
by + Yuy — 1 while decreases by 2, hence the change in Z in this case is at most the 
Oip^l'^'p^ jQg2 jj^g ^]^g error-terms in our approximation for the degrees and co-degrees. The 
probability that we choose a triangle that does not contain u yet includes an edge in T„ is 0{p'^) and 
selecting such a triangle changes Z by 0(1). The choice of any other triangle changes Z by 0{p^) 
due to the change in p. Altogether, \AZ\ = 0{n^l'^p^ log^n) and E[(AZ)2 | = Oip^"/ log^ n). 
Let po = p{io) = 1 — Gio/n^ and recall that the definition of io and our bound on \ AZ\ ensure that 

Z{io) < (/i-;u)n3/2(po)i+^<I> + 0(ni/2(po)7iog2n) < -^fin^/\po)^+nogn , 

where the factor of | absorbs the O(n^/^(po)^ log^ n)-term since npolog"^ n > n^~°^^^p* > y/n for 
large enough n. It then follows the supermartingale Sj = Z{{iQ + j) A tt) satisfies 

So < -i/in3/2(po)i+7iogn, 

max \ Sj+i — Sj\ = O (^n^^^p^ log^ , 

5^E[(5,+i - Sjf I J-j] = O (n'ipoy+'nog'n 
3 

where T'j = J^io+j- Here \So\ maxj \Sj^i — Sj\ = 0(71^(^0)^'''^^ log^ n) due to one extra log factor 
between these expressions and therefore Theorem 13.11 establishes that for some fixed c > 

p(u,{s, > 0}) < exp{-e ;;^j2Plog^!! } = {— s '-"^^ 

We conclude that w.h.p. X{i) < ^n^/'^p^'^^<^ for all u and i, and the same argument shows that 
w.h.p. X{i) > —^n^/^p^^ for all u and i. This concludes the proof of Eq. (j2.1ip . 
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3.5 Vertex degrees (1^) 

The analysis of the degrees will be straightforward using our estimate (|2.11|) for Tu, the number of 
inner edges in the neighborhood of a vertex u, since Yu changes iff the triangle selected is of the 
form uxy (in which case it decreases by 2). Indeed, setting 

X{i) = Yu-np, 

our bounds on and Q imply that 

2T 6 _ py,2^2/in3/V+7^ 6 

JH. £1A y-j — — H — — -j — o o ox 1 :t.o "I 



where the term pX"^ / {tt^p^) was absorbed into the 0(p^'^~'^<I>^/n^)-term since \X\ < K^/np'f-'^^ 
by (|2.10p . and this latter error-term was thereafter omitted as it is o{Ap^~'^^/n^^'^) by (j3.23p . 

Now consider the upper critical interval for given by 

= (/tn^/ V"^^ , Kn^/ V"^^) where k = (l - log"^ n) k , 



and as before assume that io is the first time at which X(io) > kn^^'^p^ ^<I> and define the stopping 
time Ty = min{i > ■ X{i) < kn^^'^p^^^^}. With these definitions, the variation in the variable 

Z{i) = X{i) - Kn^/V"^^ 

consists of AX as weh as A [-Kn^/V~^^] < (6 + 0{n~^)) - l)Kp^-^n-^/'^<^ + 0{n^'^/'^), hence 

E [AZ I , 4 > i] < [ - 2k + 2^ + (7 - 1)k + 0(A)] ^^3^^ < [ - (3 - 7)^ + 2fi] ^3^^ , 

where we used the fact that 7 — 7 = 6/ log n eliminates the term 2(k — k) = 2/ log n and 0(A)-term 
for large n. Hence, Sj = Z{{io + j) A Ty) is a supermartingale as long as 

(3-7)k>2^. (3.29) 

Furthermore, in each round we either select a triangle incident to u, an event which has probability 
0{l/n) and changes Z by 0(1), or we do not affect Yu and thus change Z by 0(l/n) due to the 
change in p. Thus, \A.Z\ = 0(1) while E[(AZ)^ | J^i] = 0(l/n), and we conclude that for large 
enough n the supermartingale Sj has the following attributes: 

So < -iAjni/2(po)7-i log^ ^ inax \Sj+i - S,\ = 0(1) , 5^E[(5j+i - S,)^ \ J-j] = O(npo) , 

j 

where po = p{io) = 1 — 6io/n^ and = ^i^+j- Since y/n{pof'~^' > y/n{p*)'^~'^ > clog^ n we deduce 
that |5o| maxj \Sj+i — Sj\ = o(npo) and thus Theorem 13.11 establishes that for some fixed c > 

P(U,{5, > 0}) < exp|-c "(^°^'^^[^°g'" } = {-cpr^log^n} < e-'^^^". 

Altogether, w.h.p. X{i) < nn^^^p^"^^ for all u and i, and similarly X{i) > —Kn^^^p^^^^ w.h.p. 
for all u and i. This completes the proof of Eq. (j2.10p . 
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3.6 Co-degree of triples {Yu^v,w) 

We will prove the following result from which (j2.12p will readily follow: 



Yu,v,w = np^ lb y np^ log^ n for all n, v, w and p > pi := n log^''^ n . (3.30) 



Indeed, since np^ > \J np^ log^ n for all p > pi the above would imply that Yy^^y^^ < 2np^ for all 
p > Pi- Since Yu^v,w is monotone decreasing along our process whereas 2npi = log^ n it would then 
follow that Yii^y^yj < maxjnp^ + \/ np'^ log^ n, 2npi} for all p, as required. 

It remains to prove (j3.30p . Define 

X{{) = Yu,v,w{'i) - np^- 

We have 

mx\H-- ^ + -- 



n^p^/6 \ J n 



18,, eap'^^ ^ /ap'^-^/^^log' 



2n 



T1?P V?l'^ \ V? j 

Now suppose that io is a first round at which X(io) > ^\/np^ log^ n, that is X enters the critical 
interval 



ly = ( 3 Y np^ log n , y np^ log n ) . 



Setting Ty = min{i > iq : X{i) < | \/ np^ log^ n} and 

Z(i) = — Y^^i^J^^og^ 

we get 

^r, , ... 12pi/2 log^/2 n 6ap')'^> ^/apT'-3/2$iog2^\ fog^/s „ 

E[AZ\T,,Ty>i\ < — + _^ + -2 + 



^3/2 \ 71,2 y n'^/2 



3/2 



Since 7 = 2 — \/2 — o(l) > ^ for large n, it follows that p^^ = 0{p^^'^ log^ n) = o{p^/^ log^''^ n) and 
hence the first term in the above r.h.s. is equal to —(3 + o(l))(p ^°^3"' )''^^^. As for the second term 
there, recall that p > p* and so by (j2.13p we have p^~'^^ > 6a^<&^n~^ log^ n. In particular, 

p7-3/2$lQg2^ ^plogS^y/^ $ ^ /^plog^n^ -^/^ 



7j2 I / „l/2^2-T'l — ,1/2 „ — \ n3 



^l/2p2-7 lQgl/2 7^ V n3 y log7/2 ^ ' 



and altogether we conclude that 

a/2log5/2^ 
j^3/2 



E[AZ\Ti,T^>i] <-i3 + o{l) f'~'l%"'' <0, 
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where the last inequahty holds any for sufficiently large n and confirms that Sj = Z({io + j) A Ty) 
is a supermartingale. Moreover, the equation that specified E[AX | J^i] shows that if the triangle 
selected goes through u,v or w and a vertex in Nu^v,w, which happens with probability 0{p'^/n), 
then the change in Z is 0(1), and otherwise the change in Z is 0{p'^/n). Hence, \AZ\ = 0(1) while 
E[(AZ)^ I Ti] = 0{p'^/n), thus the supermartingale Sj has the following attributes: 

So<-yn{po)nog'n , max - = 0(1) , E [(5,+i - S,)^ | J-j] = O(n(po)') , 

j 

where the factor | in the first expression (as opposed to |) treated the potential 0(1) deviation of 
Z(zo) from the lower endpoint of the critical interval. Noting that Po ^ Pi = ( ^° „ " ) ^-i^d hence 
|5o| maxj \ Sj-^-l — Sj\ = O[n{po)^), an application of Theorem 13.11 vields that for some fixed c > 0, 

P(U,{5, > 0}) < exp |-c ''^^°|2)f " I = exp (-clog^n) . 

By the usual union bound over vertices and rounds we now conclude that w.h.p. X{i) < \l np^ log^ n 
for all ti, w and i, and similarly > — \J np^ log^ n w.h.p. for all u, v, w and i, thus completing 
the proof of Eq. ([21^ . 



3.7 Number of triangles (Q) 

In the companion paper |5] it is shown (see Theorem 2 there) that w.h.p. 

Q(0 <-^ + ^n^p 

6 

throughout the process, hence it only remains to prove the lower bound in (|2.5|) . Let 
and recall that due to (12.21) we have 



xyz£Q xy£E 

To bound E[AX | J^i] from below we will thus need an upper bound on '^xyeE^x,y Recall that 
^xyeE^^,y ~ '^Q ^^^^ ^^,y ~ "-^'^ ^ an^/^p^$ by Eq. (|2.6|) . hence we can apply Lemma 
together with the fact that |-E(i)| = n'^p/2 — n/2 to obtain that 



< ^ + 0(p) + (12 + 0(A))aV^-'<J>' = ^ + (12 + 0(A))aV^-'^' , 

where in the last equality we absorbed the 0(p)-term into the 0(A) error-term factor of the last 
expression since the facts 7 < 1 and A$^ > clog^ n imply that = o(A$^), i.e. Ap^^~^<l>^ — t- 00. 

Adding this to our estimate for E[AQ | Ti] while observing that A{—^n^p^) = 2>np'^ + 0{p/n) yields 

E[AA: I F^] > - (12 + 0(A))aV^^^^^ + 3n/ + 0{p/n) . 

n'^p 
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As before we incorporate the 0{p/n) term into the 0(A) error and using the definition of X we 
can then rewrite the above as an upper bound on A(— X), as follows: 

E[A(-A) I J-,] < -^X + (12 + 0(A))aV^~^^^ • 
Now assume that is the first round where X drops below —a^nPp^'^"^^'^, i.e. enters the interval 

Iq = (^-a^nV^"^^^ , - a^nV^~^^^) where a = (l - log"^ iif^^ a . 
Further let tq = min{i > io : X{i) > — Q^n^p^'^~^$^} and 

Z(i) = -X(i) - Q2nV^"i$2 
Since A [-a^jiV"^^^^^] < (6 + 0{n^^)) {2j - l)a2p27-2^2^ ^^le upper bound on A(-A) gives 
E[AZ\Ti,TQ>i]< [- Sa^ + 2a^ + (27 - l)a'^ + 0(A)] Gp^^^^^^ < -12(1 - 7)0 V^"^^^ , 

where the last inequality used the term 7 — 7 = 6/ log n to both cancel the 0(A)-term and replace 
by 3a^. As 7 < 1 we deduce that Sj = Z{{iq + j) A tq) is indeed a supermartingale. 

Next consider the one-step variation of Z. Denoting the selected triangle in a given round by 
xyz, the change in Q following this round is at most Yx^y + Yx,z + Yy,z and in light of our co-degree 
estimate (j2.6p this expression deviates from its expected value of 3np^ by at most 3an^/^p^ log^ n. 
In particular, |AZ| = 0(-v/n(po)'^ log^ and letting po = p(^o) = 1 — Gio/n^ this ensures that 

Z{i^) < {c? - a^)n\p^f'^-^<^^ + 0[V^{poV log2 n) < -la^n'ip^f'^-^ log^ n , 

where the last inequality holds for large n since ?i'^/^(po)^~^ log n tends to 00 with n. Recalling that 
there at most n^po steps remaining until the process terminates, Hoeffding's inequality establishes 
that for some fixed c > 0, 



{Uj{Sj > 0}) < exp <^ -c ^ = exp <^ -cnpo' log n \ < e 



n 



We conclude that w.h.p. X{i) > —a'^n?p^'^ for all i, completing the proof of Eq. (12. 5p . 

This completes the proof of Theorem 12. 1[ □ 
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